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ABSTRACT 

Context. 2MASSJ05352 184-0546085 (2M0535-05) is the only known eclipsing brown dwarf (BD) binary, and so may serve as a 
benchmark for models of BD formation and evolution. However, theoretical predictions of the system's properties seem inconsistent 
with observations: /. The more massive (primary) component is observed to be cooler than the less massive (secondary) one. //. The 
secondary is more luminous (by ^ 10^^ W) than expected. Previous explanations for the temperature reversal have invoked reduced 
convective efficiency in the structure of the primary, connected to magnetic activity and to surface spots, but these explanations cannot 
account for the enhanced luminosity of the secondary. Previous studies also considered the possibility that the secondary is younger 
than the primary. 

Aims. We study the impact of tidal heating to the energy budget of both components to determine if it can account for the observed 
temperature reversal and the high luminosity of the secondary. We also compare various plausible tidal models to determine a range 
of predicted properties. 

Methods. We apply two versions of two diff'erent, well-known models for tidal interaction, respectively: /. the 'constant-phase-lag' 
model and //. the 'constant-time-lag' model and incorporate the predicted tidal heating into a model of BD structure. The four models 
diff'er in their assumptions about the rotational behavior of the bodies, the system's eccentricity and putative misalignments if/ between 
the bodies' equatorial planes and the orbital plane of the system. 

Results. The contribution of heat from tides in 2M0535-05 alone may only be large enough to account for the discrepancies between 
observation and theory in an unlikely region of the parameter space. The tidal quality factor Qbd of BDs would have to be 10^^ and 
the secondary needs a spin-orbit misalignment of > 50°. However, tidal synchronization time scales for 2M0535-05 restrict the tidal 
dissipation function to log(2BD) ^ 4.5 and rule out intense tidal heating in 2M0535-05. We provide the first constraint on Q for BDs. 
Conclusions. Tidal heating alone is unlikely to be responsible for the surprising temperature reversal within 2M0535-05. But an 
evolutionary embedment of tidal eff'ects and a coupled treatment with the structural evolution of the BDs is necessary to corroborate 
or refute this result. The heating could have slowed down the BDs' shrinking and cooling processes after the birth of the system 
^ 1 Myr ago, leading to a feedback between tidal inflation and tidal heating. Observations of old BD binaries and measurements of 
the Rossiter-McLaughlin eff'ect for 2M0535-05 can provide further constraints on Qbd- 

Key words. Celestial mechanics - (Stars:) binaries: eclipsing - Stars: evolution - Stars: individual: 2MASSJ05352 184-0546085 - 
Stars: low-mass, brown dwarfs 



1. Introduction 

2MASS J05352184-0546085 (2M0535-05) is a benchmark ob- 
ject for brown dwarf (BD) science since it offers the rare 
opportunity of independent radius and mass measurements 
on substellar objects. The obs erved values constrain evolu- 
tionary and s tructu r al m odels (D' Antona & Mazzitelli 119971: 
BaraflTe et alJ Il998l: Ichab rier & Baraffe 2000; BaraffeeLaLl 



both components formed together, as commonly believed, then 
this system allows for eflTective temperature (Tq^) and luminosity 
(L) measurements of two BDs at the same age. 

However, this system is observed t o have an unexpected tem- 
perature reversal (IStassun et al.ll2006b . contravening theoretical 
simulations: the more massive component (the primary) is the 
cooler one. From the transit light curve, the ratio of the eflTec- 
tive temperatures can be acc urately determined to Tew.2 1 Tew. i = 

L050 + 0.002 (iMohantv et al.ll2009l:[G6mez Maqueo Chew et alJ 

Nebulae, a star-forming region with an age of 1 (±0.5) Myr. If l2009l) . From spectroscopic measurements then, the absolute 



20021: Ichabrier et al.ll2007F . 2M0535-05 is located in the Orion 
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values can be constrained. The p rimary, predicted to have 
Teff,! ^ 2 870K (lBaraffeetal.1119981) . has an observed value of 
^ 2 TOOK, whereas the surface temperature of the secondary, 
predicted to be Teff 2 ~ 2 750 K, is most compatible with Teff 2 ~ 
2 890K. 

One explanation for the temperature discrepancies is sup- 
pression of convection due to spots on the surface of the pri- 
mary. If a portion of a BD's surface is covered by spots, its 
apparent temperature will be reduced, resulting in an increase 
in the estimated radius in order for the measured and expected 
luminosities to agree ( Chabrier et al. 2007). With a spot cov- 
erage of 30 - 50% and a mixing length parameter a = I 
most of the mismatches between predicted and observed radi i 
for low-mass stars (LMS) can be explained (iRibas et al.ll2008l) . 
Observations of spots on both of the 2M0535-05 components 
(lG6mez Maqueo Chew et al. 2009), as inferred from periodic 
variations in the light curve, and measurements on the Ha 
line of the combined spe ctrum during the radial velocity max- 
ima (iReiners et al ] l2007l) suggest that enhanced magnetic ac- 
tivity and the accompanying spots on the primary indeed play 
a key role for its temperature deviation. But even if the spot 
coverage on the primary serves as an explanation for the pri- 
mary's reduced Teff, the seconda ry's luminosity over shoot of 
^ 2.3 • lO^^^ W, as compared to the lBaraffe eTaP (Il998h models, 
suggests some additional processes may be at work. 

The temperature reversal between the primary and secondary 
may result from a difference between their ages. The secondary 
could be ^ 0.5 Myr older than the primary, as proposed by 
IStassun et all ^07) (see also D'Antona & Mazzitelh (1997)). 
A difference of 0.5 Myr could allow the secondary to have 
converted the necessary amount of gravitational energy into 
heafl, which would explain its luminosity excess. But evolution- 
ary n i odels a re very unc ertain for ages ^ 1 Myr (Baraffe et al.' 
I2002t IWuchte rl 2005: Marlev et al.ll2007l: [Mohanty elal, 2002) 
and, in any case, the age determination and ph ysical natures 
of these ver y young objects is subject to debate (IStassun et al.l 
I2008ll2009h . Furthermore, the mutual capture of BDs and LMS 
into binary systems after each component formed indepen- 
dently is probably too infrequent to account for the large num- 
ber of eclipsing LMS binaries with either temperature reversals 
or inflated r adii (Guenther et al.""200l'; 'Coushhn & Shaw'(2007|; 
iRibas et al.i r2008: Cakirh et al. 2009; Morales et al. 2009). 

Here, we consider the role that tidal heating may play in de- 
termining the temperatures of the BDs. In Table [T] we show the 
parameters of 2M0535-05 necessary for our calculations. The 
computed energy rates will add to the luminosity of the BDs in 
some way (Sect. 12.31) and will contribute to a temperature devi- 
ation compared to the case without a perturbing body (Sect. [3]). 
All these energy rates must be seen in the context of the lumi- 
nosities of the BDs: Li ^ 8.9-10^'^ W (luminosity of the primary) 
and L2 ~ 6.6 • 10^"^ W (luminosity of the secondary). At a dis- 
tance a to the primary component, its luminosity is distributed 
onto a sphere with area 47ra^. The secondary has an effective - 
i.e. a 2D-projected - area of tiRI. With Fi^^ as the flux of the 
primary at distance a, the irradiation from the primary onto the 
secondary Li^2 is thus given by 

L 

Li^2 =7:Rl Fi,a = 7:Rl =^17^. (D 



^ In contrast to the Baraffe et"aLl (Il998h tracks, the models by 
iD'Antona & Mazzitellil ( 1997h predict a temperature increase in BDs 
for the first ^ 30 Myr of their existence. 



Using that equation, we calculate the mutual irradiation of the 
BDs: Li^2 ~ 8.5 • lO^^ W and L2^i ^ 1.0 • lO^^W. These en- 
ergy rates are two and three orders of magnitude lower, respec- 
tively, than the observed luminosity discrepancy. Hence, we as- 
sume that mutual irradiation can be ignored. This simplification 
is in contrast to the cases of the potentially inflated transiting 
extrasolar planets WASP-4b WASP-6b WASP - 12b, and TrES- 
4, where stellar irradiation (llbgui et al. 1 12009^ dominates tidal 
heating by several magnitudes. 

Various tidal mod els haven been used to calculate tidal 
heating i n exoplanets (jBodenheimer et al.l 1200 H: iJackson et al.l 
!2008aBbl: I Barnes et al. 1 12009"), which may in fact be responsible 
for previous discrepan cies between interior m odels and radii of 
transi ting exoplanets (iJackson et al.l l2008allbl : llbgui & Burrowsl 
I2OO9I) . This success in exoplanets motivates our investigation 
into BDs. While many different tidal models are available, there 
is no consensus as to which is the best. For this reason, we 
apply a potpourri of well-established models to the case of 
2M0535-05 in order to compare the different results. As we 
show, tidal heating may account for the temperature reversal and 
it may have a profound eflTect on the longer-term thermal evolu- 
tion of the system. 

The coincidence of Porh/Pi ~ 2.9698 ^ 3, with Porb as the 
orbital and Pi primary's rotation period, has been noted before 
but we assume no resonance between the primary's rotation and 
the orbit for our calculations. These resonances typically occur 
in systems with rigid bodies where a fixed deformation of at least 
one body persists, such as in the Sun-Mercury configuration with 
Mercury trapped in a 3/2 spin-orbit resonance. We assume that, 
in the context of tides, BDs may rather be treated as fluids and 
the shape of the body is not fixed. 

With this paper, we present the first investigation of tidal in- 
teraction between BDs. In Sect. [2] we introduce four models for 
tidal interaction and discuss how we convert the computed en- 
ergy rates into an increase in eflTective temperature. Sect. [3] is 
devoted to the results of our calculations, while we deal with the 
observational implications in Sect. IH We end with conclusions 
about tidal heating in 2M0535-05, and in BDs in general, in 
Sect.E 

2. Tidal Models 

Two qualitatively diflTerent models of tidal dissipation and 
evolution have been developed over the la s t century: The 
^constant-phase-lag' (iGoldreich & Sotej 1 19661: IWisdomI l2008l: 
lFerraz-Melloetal.l l2008l Wis08 and F M08 in the following), 
and the 'constant-time-lag' model (lHutlll98ll Hut81 in the fol- 
lowing). In the former model, the forces acting on the deformed 
body are described by a superposition of a static equilibrium po- 
tential and a disturbing potential (FM08). The latter model as- 
sumes the time between the passage of the perturbing body over- 
head and the passage of the tidal bulge is constant. Although 
both models have been used extensively, it is not clear which 
model provides a more accurate description of the effects of 
tides, so we apply formulations of both models. 

In the 'constant-phase-lag' model of FM08, quantitative ex- 
pressions have been developed to second order in eccentric- 
ity e while the others include also higher orders. Higher and 
higher order expansions require assumptions about the depen- 
dence of a body's tidal response to an increasing number of tidal 
frequencies, which involves considerable uncertainty. Therefore 
higher order ex pansions do not necessarily provide more ac- 
curacy (FM08: iGreenbergl l2009l) . In the 'constant-phase-lag' 
model of Wis08, expressions in e are developed to 8* order. 
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The 'constant- time-lag' model of HutSl does not include pos- 
sible obliquit i es, wh ile an enhanced version of that model by 
iLevrard et all (l20Q7h (Lev07) does. 

Tidal dissipation in BDs has not been observed or even con- 
sidered previously, and hence, neither model should take prece- 
dence when calculating their tidal dissipation, especially since 
neither tidal model is definitive ( Greenberg 2009). As our in- 
vestigation is the first to consider tidal eff'ects on BDs, we 
will employ several applicable, previously published models to 
2M0535-05. By surveying a range of plausible models and in- 
ternal propertie s, usually encapsulated in the 'tidal dissipation 
function' Q ( G oldreich & Soter 1966), we may actually be able 
to determine which model is more applicable to the case of BDs 
- assuming, of course, that tidal dissipation contributes crucially 
to the observed temperature inversion. 



Table 1. Orbital and physical parameters of 2M0535-05 



PROPERTY 


OBSERVED VALUE 


a, semi-major axis^ 


0.0407 +0.0008 AU 


e, eccentricity^ 


0.3216 +0.0019 


Porb, orbital period^ 


9.779556 +0.000019 d 


/, orbital inclination to the line of sight ^ 


88.49 +0.06° 


age^ 


1 +0.5 Myr 


Teff,!, primary effective temperature^ 


2715 +100K 


7'eff,2/7'efF,i, effective temperature ratio ^ 


1.050 +0.002 


Ml, primary mass^ 


0.0572 +0.0033 


M2, secondary mass^ 


0.0366 +0.0022 


Ri, primary radius^ 


0.690 +0.01 IPo 


R2, secondary radius^ 


0.540 +0.009 Po 


Li , primary luminosity^ 


8.9-1024 +3- 10^4 W 


L2, secondary luminosity^ 


6.6-1024 +2- 10^4 W 


Pi, rotational period of the primary^ 


3.293 +0.001 d 


P2, rotational period of the secondary^ 


14.05 +0.05 d 


f efF,i , modeled T^s for the primary^ 


2 850K 


^efF,2, modeled T^s for the secondary^ 


2 700K 


Ri , modeled radius for the primary^ 


0.626 Po 


R2, modeled radius for the secondary^ 


0.44 Po 



^ Gomez Maaueo Chew et al. (2009) . ^ Baraffe et al. (1998) . 
^ assuming an uncertainty of 200 K in Teff,! and reff,2 



2.1. Constant phase lag 
2.1.1. Tidal model #1 

The potential of the perturbed body can be treated as the super- 
position of periodic contributions of tidal frequencies at differ- 
ent phase lags and the expression for the potential can be ex- 
panded to first order in those lags (FM08). Those phase lags 
^k,i I yi=o, 1,2,5,8,9 of the body that we will need for our equations 
are given by 



Qi £o,i 


= i:(2Q/ 


- 2n) 


Qi sij 


= i:(2Q,- 


-3n) 


Qi S2,i 


= i:(2^i 


-n) 


Qi S5,i 


= ^(n) 




Qi ^8,/ 


= ^(^i - 


2n) 


Qi £94 


= 





/e{l,2}. 



(2) 



where is the algebraic sign of x, thus = -h 1 V - 1 , 
n = 27r/Porb is the orbital frequency and = 27i/Pi are the 
rotational frequencies of the primary (/ = 1) and secondary (/ = 
2), Pi being their rotational periods. The tidal frequencies are 
functions of the tidal quality factor Q of the deformed object, 
which parametrizes the object's tidal response to the perturber. 
It is defined as 



— r 

27lEo Jo 



dt 



dE 
'~dt 



(3) 



where Eq is the maximum energy stored in the tidal distor- 
tion and the integral over the energy dissi pation rate -dE/dt 
is the energy los t during one o r bital c ycle (iGoldreich & Sotej 
Il966h . Although [Ogilvie & Lid (l2004l) conclude that tidal dis- 
sipation rates of giant planets are not adequately represented 
by a constant Q-value, many parameterized tidal models rely 
on this quantity. Measurements of the heat flux from Jupiter's 
moon lo during the fly-by of the Voyager 1 spacecraft, combined 
with a specific model of the history of the orbital resonance, al- 
lowed for an estimate for the quality factor Qq^ of Jupiter to be 
2_10f < < 2 • 10^ (lYodedll979l) while lAksnes & FranklinI 
(l200lh used historical changes in lo's orbit to infer that Qq^ 
is around 10^^. However. iGreenberg et al.l (l2008l) pointed out 
that Q = 00 is not ruled o ut (see also lPeale & Greenberpll98Ql: 
lloannou & LindzenI 19931) . Tides raised by Neptune on its moons 
help to constrain the planet's quality factor to 10^^^ < 2g < 
104-56 (IZhang & Hamilton! l2008h . For M dwarfs, Qm is as- 
sumed to be of order 10^, whereas for rigid bodies li ke Earth 
20 < 2 < 500 (iRav et al.ll200ll: iMardling & Linll2004 and ref- 
erences therein). For BDs, however, Q is even more uncertain, 
thus we will handle it as a free parameter in our procedures. 

FM08 allows for the tidal amplitude to be different from 
what it would be if the tide-raising body were fixed in space. 
This concern is met by the dynamical Love number under the 
assumption that the tidally disturbed body had infinite time to re- 
spond. Without better knowledge of a body's response to tides, 
we assume the dynamical Love number is the same as the poten- 
tial Love number of degree 2, ^2- For the g as planets of the solar 
system , this number has been calculated by lGavrilov & Zharkovl 
(1 19771) . BDs may rather be treated as poly tropes of order n = 3/2 
(L Baraffe, private communi cation). We infer the L ove number 
from the relation k2 = 2^aps (iMardling & Lin|[2002b and use the 
tables of apsidal motion constants ^aps given in lBrooker & Olid 
(1955). These authors provide numerical calculations for ^aps for 
a polytrope of n = 3/2. We find k^p^ = 0.143 and thus = ^2 = 
0.286. This places for BDs well in the regime spanned by 
the gas giants of the solar system: Jupiter (k2 = 0.379), Saturn 
(k2 = 0.341), Uranus (fe = 0.104) and Neptune (fe = 0.127) 
( Ga vrilov & Z harkov 1977). 

Before we proceed to the equations for the tidal heating rates, 
we sum up those for the orbital evolution of the semi-major axis 



4 



Heller et aL: Tidal effects on brown dwarfs: Application to 2M0535-05 



a, the eccentricity e and the putative obliquity i//. The latter pa- 
rameter is the angle between the equatorial plane of on e of the 
two b odies in a binary system and the orbital plane ( Win n et aP 
l2005l) . frequently referred to as spin-orbit misalignment. We use 
Eqs. (56), (60) and (61) from FM08 but our equations for a bi- 
nary system with comparable masses need slight modifications 
since both constituents contribute significantly to the evolution 
of a and e. We add both the terms for the secondary being the 
perturber of the primary (/ = I, j = 2) and vice versa, since only 
spin-orbit coupling is relevant, whereas spin- spin interaction can 
be neglected. This results in 



da ^ 3ka,iMjR^^n 



dt 



147 1 



-3s5,i]-4S^[so,i-ss,i]\ (4) 



de ^ 
dt^~^ 



1,2 



SekdjMjR^^n 



49 1 

2soj - —Si^i + -£2,i + 3£5,i I , (5) 



d^ _ 3kajMjR^^n 
dt AMiU^ 



S i {-£Q,i + ^8,/ + -£9,i) , (6) 



where ^d,/ is the dynamical Love number, M/ the mass and Ri the 
radius of the deformed BD, S / := sin(0^/), with ij/i as the obliquity 
of the perturbed body, and \ yt=o, 1,2,5,8,9 are the tidal phase lags, 
given in Eq. Q. 

The total energy that is dissipated within the perturbed body, 
its tidal energy rate, can be determined by summing the work 
done by tidal torques (Eqs. (48) and (49) in FM08). The change 
in orbital energy of the body due to the body is given by 



Sk^iGM^.R^ 147 1 

^orb,/ = ^ ^ 4£:o,/ + e\-20soj -h —sij -h -6:2,/ 



-3^5,/] -45? [^0,-^8,/]) (7) 
and the change in rotational energy is deduced to be 



3ka,iGMjR^ 
8^ 



- _ 



( 4so,i + e^[-20so,i + 49£:i,/ -h S2,i] 
-^2S^[-2so,i-^ss,i-^S9,i]), (8) 



where G is Newton's gravitational constant. The total energy re- 
leased inside the body then is 



^#1 



orb,/ 



+ 4*ou)>0. 



(9) 



The greater-than sign in this equation is true, since either Q/ < n 
and orbital energy is converted into rotational energy, or Q/ > n 
and the body is decelerated by a transfer of rotational energy into 
orbital energy. In both cases, the dynamical energy of the system 
is released within the distorted body. For = 0, e.g., Eqs. d?]) 
and m yield E^J^ . = -p - (4 ^ 51 + AS})IQi and . = 0. 

The approach for the calculation of tidal energy rates with 
tidal model #1 depends on processes due to non- synchronous 
rotation via Sk^i = Skji^i, n) and includes a putative obliquity ij/i 
and terms of e up to the second order. After inserting the orbital 



and rotational periods for 2M0535-05, these equations reduce 
to 



^tid,l 



^tid,2 - 



— ■ 1-^ ( [4 -h 30^2]Di - [4 -h 5\e^]n ) , 

?>k^^2GM\Rl 



( [4 -h 56^2]^ + [2Sl-A- 28^2]Q2 ) • 

(10) 



Interestingly, for these particular values of Di, Q2 and n, the 5" 1- 
terms for E"*^ ^ cancel each other, so that it is not a function of 
j/^i, whereas E"*^ ^ does depend on \l/2. 



2.1.2. Tidal model #2 

The model of Wis08 includes terms in eccentricity up to the 8* 
order, predicting higher tidal energy rates than for the equations 
of model #1. Equations for the evolution of the orbital parame- 
ters are not given in Wis08. Furthermore, in his theory the per- 
turbed body is assumed to be synchronously rotating with the 
orbital period. Since this is not the case for either of the BDs 
in 2M0535-05, the following equations will only yield lower 
limits for the tidal heating. The tidal heating rates are given by 



^tid,/ 



2lk2,iGM]R^.n 



2Qia^ 



(11) 



with 



o fHut A fHut 1 f Hut 

^w.s(e,if^i) - --^ Ci + --^ (1 + C,) 

3 ^2^Wis 

where we used C, := cos(^,) and 



(12) 



p = 






/•Hut 
Jl - 




255 4 185 


/•Hut 
Jl - 






/•Hut 
^5 - 


1 3^H 




/•Wis 
J?, - 




24 1 6 



(13) 



following the nomenclature of |Hu3 (Il98lh and lWisdonJ (l2008h 
as indicated. Furthermore, fe,/ is the potential Love number of 
degree 2 for the component of the binary system and A/ is 
a measure of the longitude of the node of the body's equator 
on the orbit plane with respect to the pericenter of its orbit. In 
order to estimate the impact of A/ in the last term in Eq. ([T2b . 
we assume this impact to be as large as possible. A/ = 0, and 
compare it to the preceding terms. We find that for the case of 
2M0535-05 the first three terms are of order 1, whereas the term 
connected to A/ varies between 10"^ and 10"^, depending on ij/i. 
These irrelevant contributions give us a justification to neglect 
the unknown values of A/ in 2M0535-05 for our computations, 
facilitating the comparisons to the other models. 
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2.2. Constant time lag 
2.2.1. Tidal model #3 

Instead of assuming phase lags and superposition of frea i iency- 
dependent potentials, the 'equilibrium tide' model by iHutI (1 1 9 8 ll) 
invokes a constant time lag r between the line joining the cen- 
ters of the two bodies and the culmination of the tidal bulge on 
the distorted object. With that assumption, the model of HutSl 
is mutually exclusive wit h the assumption of a fixed angle lag 
(iGoldreich & Soterl 19661) : in general, a fixed time lag and a fixed 
angle lag result in very diff'erent behaviors of the tidal bulg43- As 
for the case of the 'constant- time-lag' model, we first sum up the 
equations governing the behavior of the orbital evolution. With 
the purpose of easing a comparison between HutSl 's equations 
(Eqs. (9), (10), and (11) therein) and Eqs. © - Q from this 
paper for the theory of the 'constant-phase-lag' model #1, we 
transform the former into 



where 



da ^ 
d7 ~ ^ 

i=l,2 



-6k,^,jGMjR^ 



Mj 
' Mi 



de_y -21k^^,^iGMjR\e 
dt ~ .^^ 

I - 1,2 



2/3^ 



1 -h — ^ 

Mi 



pl5 



fHut 
h 



n 



(14) 



18 n 

(15) 



jsi2 a- 



M, 



\ a I n 



(16) 



■Hut 



pl5 



fHut 



(19) 



Unfortunately, with these equations for the tidal energy rates 
model #3 neglects a potential obliquity of the body, which pre- 
vents us from a direct comparison with the other tidal models. 

2.2.2. Tidal model #4 

Lev07 extended HutSl 's formula for the tidal energy rate to 
the case of an object in equili brium rotatior0 and they include d 
possible obliquities (see also iNeron de Surgy & Laska^ll997l) . 
though they do not give the equations for the orbital evolution. 
LevOT's equations are equivalent to 



^tid,/ ~ 



3k2,iGMjR^n 



(20) 



where 



pl5 



fHut/ ^9 



1 -h 



l-2/Sl 



(21) 



The 'annual tidal quality factor' is given as 2;^^ = n r. Even 
though LevOT's equations invoke Qn and their equations resem- 
ble those of the models with constant phase lag, their approach 
still assumes a constant- time-lag. Since Lev07 do not explicitly 
connect their to the Q of FMOS (model #1) and WisOS (model 
#2), we keep Q and Qn as two diff'erent constants for our further 
treatment. 

With these expansions, Eq. (l2Ql) involves terms in eccentric- 
ity up to order e^. But since model #4 assumes tidal locking, i.e. 
E^^ is not a function of this model a lso yields just a lower 
limit for the heating rates (IWisdomll200"8h . 



with /^aps,/ as the apsidal motion constant of the perturbed body 
(see Sect. 12.1.11) , rl . as the radius of gyration of the body, 
which is defined by the body's moment of inertia // = M^r^ .R^, 
and 



/Hut = 1^11,2^ 4 ^ 5 6 

4 8 64 



/Hut ^ 1 + + 



1 



(17) 



2.3. Converting tidal heating into temperature increase 

Now that we have set up four distinct models for the calcula- 
tions of the additional tidal heating term for the BDs, there are 
two physical processes that will be driven by these energy rates: 
tidal inflation and temperature increase. Let's take L as the lumi- 
nosity of either of the two 2M0535-05 BDs that it would have 
if it were a single BD and R and f etf as its corresponding radius 
and eff'ective temperature. The n, by the Stefan-Boltzmann law 
(IStefanll 1 8791: iBoltzmannll 1 881 



Hut8 1 then calculates the energy dissipation rate within a binary 
system, caused by the influence of one of the two bodies on the 
other, as the change in the total energy E = £^orb + ^^rot- Here, 
£'orb and E^ot are the orbital and rotational energies of the body 
(Eqs. (A28) - (A35) in Hut81). For the tidal heating rates of the 
constituent within the binary, this yields 



3h^,jGMjRy 



^#3 
"tid,/ 



(18) 



^ If ^ = and if/ - 0, then there is a single tidal lag angle s and 
the tidal dissipation funtion can be written as 2 = \/s = l/irn). For 
the course of an orbit, where the tidal evolution of n is negligible, both 
Q and r can be fixed. However, in a general case where r is constant 
in time, Q will decrease as the orbital semi-major axis decays and n 
increases. So Q would not be constant. 



L = 47r^VsBf; 



eff' 



(22) 



where osb is the Stefan-Boltzmann constant. The radial expan- 
sion in the binary case is given by dR = R-R and the temperature 
increase by dT = T^^ - fetf • In its present state, the BD has a 
luminosity 



(23) 



where Ei^ is some additional internal energy rate. Solving Eq. 
I for the temperature increase yields: 



dT 



I 4;r7? VsB 



'eff 



1/4 



■ ^eff • 



(24) 



WisOS calls this 'asymptotic nonsynchronous rotation'. 
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In the next step, we quantify the amount of tidal energy that 
is converted into internal energy, leading to an increase in ef- 
fective temperature. Since we will use the virial theorem for an 
ideal, monoatomic gas to estimate the partition between inter- 
nal and gravitational energy, we first have to assess the adequacy 
of treating the 2M0535-05 BDs as ideal gases. We therefore 
show the degeneracy parameter ^ = k^T l(k^TY) as a function 
of radius in Fig.[T](Chabrier & Baraflte 2000, 1. Baraff'e, private 
communication.). Here, is the Boltzmann constant, T is the 
local temperature within the gas and = ^b^f is the Fermi en- 
ergy of a partially degenerate electron gas with an electron Fermi 
temperature Tp. With respect to M, Teff and log(g), g being the 
body's gravitational acceleration at the surface, the BD struc- 
ture model corresponds to that of the primary, but with an age of 
4.9 Myr. We find that for most of the BD, i.e. that portion of the 
structure in which the majority of the luminosity is released, ^ 
is of order 1 . This means that we may indeed approximate the 
BDs as ideal gases. 

With the time de rivative of the virial theore m for an 
ideal monoatomic gas (iKippenhahn & Weiger 3 fT99Ql Sect. 3.1 
therein). 



L - Ei^ - -Eq/2, 



(25) 



where £"0 is the temporal change in gravitational energy, we find 
that half of the additional tidal energy is converted into internal 
energy and the other half causes an expansion of the BD. There 
are currently no models for tidal inflation in BDs and the treat- 
ment is beyond the scope of this paper. Instead of including the 
modeled BD radii Ri into Eq. ([24b we avoid further uncertainties 
and fixR/R = 1 (see Sect. [3 for a discussion of tidal inflation in 
the evolutionary context). The increase in efl'ective temperature 
due to tidal heating then becomes 



dT 



I Em/2 
\47r7?VsB 



1/4 



f eff- 



(26) 



For Tq^j we took the values predicted by the lBaraff'e et al.l (1 19981) 
models (see Table [T]). 

Our neglect of tidal inflation makes this temperature in- 
crease an upper limit. Given that this neglect is arbitrary, we 
estimate how our constraints for log(22) = 3.5 and 1/^2 = 50° 
would change, if tidal inflation played a role in 2M0535-05. 
Comparing the observed radii of both BDs with the model pre- 
dictions (see Table [B, radial expansions of 10% for the pri- 
mary and 20% for the secondary seem realistic. Theoretical 
investigations o f tidal heating on the in flated transiting planet 
HD209458b by (llbgui & BuiTowsll2009l) support an estimate of 
tidal inflation by 20%. As a test, we assumed that the sec- 
ondary BD in 2M0535-05 is tidally inflated, where its radius 
in an isolated scenario would be 80% of its current value, i.e. 
R = 0.8 R in Eq. (|24l). In the non-inflated scenario with R/R = 1, 
the BD would reach a temperature increase of dr = 60 K at 
log(e2) = 3.5 and 1/^2 = 50° with model #2 (see Sect. [33]). With 
the inflation, however, log(22) ~ 2.7 is needed to achieve the 
same heating at 1/^2 = 50°, whereas no obliquity at log(22) = 3.5 
would yield significant heating. Thus, if tidal inflation in the sec- 
ondary BD increases its radius by 20%, then the value for the dis- 
sipation function required to yield the same Tq^ would be about 
0.8 smaller in log(2) than in the case of no inflation. Therefore, 
the temperature we report in Sect. [3]may, at worst, correspond to 
log(G) that is smaller by 0.8. 



3. Results 

3.1. Orbital evolution 

In order to get a rough impression of how far the orbital config- 
uration of the system has evolved, we used the equations given 
in FM08, to compute the change of its eccentricity e and of a 
possible obliquity i//2 of the secondary within the last 1.5 Myr. 
Since this time span is the upper bound for the system's age, 
confined by its localization within the Orion Nebulae and in- 
dicated by comparison with BD evolutionary tracks, we thus 
get the strongest changes in e and i//2. If any initial obliquity 
would be washed out already, could be neglected in the cal- 
culations of tidal heating. Furthermore, the measured eccentric- 
ity e could give a constraint to the tidal dissipation function Q. 
Computations based on the theory of 'constant- time-lag' yield 
qualitatively similar results. 

For the evolution of e, we relied on Eq. ([5]). We took the ob- 
served eccentricity ^ = 0.3216 as a starting value and evolved 
it backwards in time. To evolve the system into the past, we 
changed the sign of the right side of the equation. Furthermore, 
we assumed that the quality factors Qi and Q2 of the pri- 
mary and secondary are equal, leading to 2i = ^2 =• G and 
^k,i •= ^k,i(^i, 2), because we are merely interested in a tenta- 
tive estimate so far. This assumption should be a good approxi- 
mation due to the similarity of the both components in terms of 
composition, temperature, mass, and radius. 

The observed eccentricity of the system might give a con- 
straint to the possible values for Q since de/dt depends on Q via 
Sk^i. Certain Q regimes could be incompatible with the observed 
eccentricity of the system at a maximum age of 1.5 Myr, if these 
Q values would have caused the eccentricity to decay rapidly to 
within this time. However, our simulations (Fig. [2]) show that 
the system has not yet evolved very far for the whole range of 
Q and that the eccentricity of 2M0535-05 is in fact increasing 
nowadays. In this system, circularization does not occur. The ob- 
served eccentricity of 0.3216 consequently does not constrain Q. 
In this first estimate, we fixed all other parameters in time, i.e. 
we neglected an evolution of the semi-major axis a, of possible 
obliquities and we used constant radii Ri and rotational fre- 
quencies We did this because we cannot yet incorporate the 
evolutionary behavior of the components' radii Ri in the context 
of tides and furthermore, there is no knowledge about possible 
misalignments ij/i between the orbital plane and the equatorial 
planes of the primary and secondary, respectively. A consistent 
evolution of Ri, however, is necessary to evolve dal&t as a func- 
tion of ij/i and j/^2, as given by Eq. ([4]). Such a calculation was 
beyond the scope of this study. 

The relative spin-geometry of the two BD rotational axes 
with respect to the orbital plane and with respect to each other 
is unknown in 2M0535-05. Anyhow, we can estimate if a pos- 
sible obliquity that once existed for one of the BDs would still 
exist at an age of 1.5 Myr or if it would have been washed out 
up to the present. We used a numerical integration of Eq. ([6]) to 
evolve j/^2 forward in time (Fig. [2]). For the secondary's initial 
obliquity 0^ini,2, we plot the state of 1I/2 as a function of the qual- 
ity factor Q2 after an evolution of 1.5 Myr. We see that even for 
a very small quality factor of 10^ and high initial obliquities the 
secondary is basically in its natal configuration today. Thus, it is 
reasonable to include a putative misalignment of the secondary 
with respect to the orbital plane in our considerations. As shown 
below, this is crucial for the calculations of the tidal heating and 
the temperature reversal. 
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3.2. Tidal heating in 2M0535-05 with model #1 

In Fig.[3l we show the results for the tidal heating rates as com- 
puted after tidal model #1 . As given by Eq. ([TOl) , the tidal heating 
of the primary does not depend on a putative obliquity, whereas 
that of the secondary does. Using this model, we find that the 
luminosity gain of the secondary is, over the whole Q range, 
smaller than that of the primary, which mainly results from the 
relation E^^^ . oc R^.. Figure [3] also shows that a growing obliq- 
uity shifts the gain in thermal energy towards higher values for a 
fixed Q2. The observed overshoot of 10^^ W in the secondary's 
luminosity can be reproduced with very small quality factors of 
22 ~ 10^ and high obliquities up to </^2 ~ 90°. 

In Fig. 131 we show the results for the temperature increase 
as per Eq. ([26b with the tidal energy rates coming from model 
#1. These rates yield only a slight temperature increase for both 
constituents. Even for low Q values of order 10"^ and high obliq- 
uities of the secondary, the heating only reaches values < lOK. 
We also see that the heating for the primary is computed to be 
greater than that for the secondary and no temperature reversal 
would be expected. If both BDs have the same Q values, then 
model #1 is unable to explain the temperature reversal. We can- 
not rule out a system in which, e.g., Qi = 10^ and Q2 = 10^, 
for which model #1 could explain the reversal. However, there is 
no reason to expect that similar bodies have Q values that span 
orders of magnitude. Hence, we conclude that model #1 can nei- 
ther reproduce the luminosity overshoot of the secondary nor the 
system's temperature reversal. 

3.3. Tidal heating in 2M0535-05 with model #2 

This model yields the highest heating rates and hence temper- 
ature increases. The contrast between the absolute energy rates 
within the primary E^^^ ^ and the secondary ^ is very small. 
In fact, for any given point in iji-Q space, the heating rates dif- 
fer only by \og{Efj^) - \og{Efjy^) ^ 0.1 (Fig. O. The 
tidal energy rates of the secondary become comparable to the 
observed luminosity overshoot at log(22) ~ 3.5 and </^2 ~ 50°, 
where E^^^^ ^ 10^^ W. A comparison of the heating rates from 
model #2 with those of model #1 for either of the BDs shows 
that model #2 provides higher rates, with growing contrast for 
increasing obliquities. 

The temperature increase arising from the comparable heat- 
ing rates is inverted for a given spot on the i//-log(Q) plane. If 
both BDs had the same obliquity and the same dissipation factor, 
the secondary would experience a higher temperature increase. 
As presented in Fig.[6l the temperature increase after model #2 is 
significant only in the regime of very low Q and high obliquities. 
Neglecting any orbital or thermal evolution of the system, the 
observed temperature reversal could be reproduced by assuming 
an obliquity for the secondary while the primary's rotation axis 
is nearly aligned with the normal of the orbital plane. We note 
that the real heating will probably be greater since model #2 as- 
sumes synchronous rotation, which is not the case for both BDs 
in 2M0535-05 (see Table[T]). The values of Q2 and i//2 necessary 
to account for the observed increase in L2 and reff,2 may thus be 
further shifted towards more reasonable numbers, i.e. Q2 might 
also be higher than 10^^ and the obliquity might be smaller than 
50°. Thus, for a narrow region in the i//-log(Q) plane, model 
#2 yields tidal energy rates for the secondary comparable to its 
observed luminosity overshoot and in this region the computed 
temperature increase can explain the observed temperature re- 
versal. 



3.4. Tidal heating in 2M0535-05 with model #3 

Since the only free parameter in this model is the putative fixed 
time lag r, we show the tidal heating rates for both the pri- 
mary and the secondary only as a function of r in Fig. [7] with 
s < T < 300 s. For this range, model #3 yields energy rates and 
temperature rises that are compatible with the observed luminos- 
ity and temperature overshoot of the secondary. For r > 100 s 
the heating rate for the secondary becomes comparable to the 
observed one, namely ^^^2 ~ 10^^ W. However, assuming a 
similar time lag t\ for the primary, the luminosity gain of the 
primary BD would be significantly higher than that of the sec- 
ondary, which is not compatible with the observations. The as- 
sumption of Ti ^ T2 should be valid since both BDs are very 
similar in their structural properties, such as mass, composition, 
temperature, and radius. 

The corresponding temperature increase is plotted in Fig.O 
It shows that the more massive BD would experience a higher 
temperature increase than its companion, assuming similar time 
lags. Since tidal heating is underway in 2M0535-05 and was 
probably similar in the past (see Sect. 12.1.11) , tidal heating af- 
ter model #3 would have been more important on the primary, 
forcing it to be even hotter than it would be without the per- 
turbations of the secondary. The temperature diff'erence between 
the primary and the secondary, which is anticipated by BD evo- 
lutionary models, would be even larger. Thus, the temperature 
inversion cannot be explained by tidal model #3. 

3.5. Tidal heating in 2M0535-05 with model #4 

The calculations based on model #4 yield significant heating 
rates in both BDs. Like in the case of models #1 and #2, the 
luminosity gain of the secondary at a fixed obliquity is, over the 
whole 2n range, smaller than that of the primary (Fig. O. As 
for model #2, the diff'erence between £*^^(</r) and is 
less pronounced than in model #1. Assuming spin-orbit align- 
ment for the primary and a pronounced obliquity of the sec- 
ondary, tidal heating rates of £*^2 - 10^^ W can be reached 
with log(2n,2) ~ 3.5 and 1^2 ~ 50°.' 

Like model #2, #4 produces a reversal in temperature in- 
crease by means of the modified Stefan-Boltzmann relation in 
Eq. ([26b , due to the comparable heating rates of both BDs and 
the significantly smaller radius of the secondary (Fig. [T0|). We 
find a reversal in tidal heating, i.e. dr2 > d^i for any given 
point in il/-Qn space. A temperature increase of > 40 K can be 
reached with log(2n,2) ~ 3.5 and </^2 ~ 50°. Since the equations 
of model #4 provide merely a lower limit due to the assump- 
tion of asymptotic non- synchronous rotation, 2n,2 might also be 
higher than 10^^ and the obliquity might be smaller than 50°. 
Similar to model #2, tidal model #4 can reproduce the observed 
temperature reversal in a narrow region of the j/^-log(2) parame- 
ter space. 

4. Discussion 

We employed several tidal models to explore the tidal heating 
in 2M0535-05. We found that, assuming similar tidal quality 
factors Q and obliquities ij/ for both BDs, the constant-phase-lag 
model #2 and the constant-time-lag model #4 yield a stronger 
increase in eff'ective temperature on the secondary mass BD than 
on the primary. For certain regimes of Q2 and j/^, the tidal energy 
rates in the secondary are of the correct amount to explain the 
larger temperature in the smaller BD. A comparison between our 
computations based on the models #1 and #2 on the one hand and 
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#3 and #4 on the other hand is difficult. The reference to a fixed 
tidal time lag might only be reconciled with the assumption of 
= ^ T as done by Lev07, which is at least questionable since 
the assumption of a fixed time lag is not compatible with a fixed 
phase lag. Furthermore, model #3 does not invoke obliquities, 
which also complicates direct comparisons of the model output. 

4.1. Constraints on the tidal dissipation function for BDs, 2bd 

4.1 .1 . The Rossiter-McLaughlin effect in 2M0535-05 

The geometric implication of the most promising tidal models #2 
and #4 is that the obliquity of the 2M0535-05 primary is negli- 
gible and that of the secondary is 1/^2 ~ 50° - provided tidal heat- 
ing accounts for the Teff reversal and the luminosity excess of 
the secondary. There does exist an observational method to mea- 
sure the geometric configuration of eclipsing systems, called the 
Rossiter-McLaughlin effect (RME) (Rossiter 1924; McLaughHn 
11924 . 

The RME appears during transits in front of rotating stars. 
Hiding a fraction of the star's surface results in the absence 
of some corresponding rotational velocity contribution to the 
broadening of the stellar lines. Thus, the changes in the line pro- 
files become asymmetric (except for the midpoint of the transit) 
and the center of a certain stellar line is shifted during a transit, 
which induces a change of the star's radial velocity. The shape of 
the resulting radial velocity curve depends on the effective area 
covered by the transiting object and its projected path over the 
stellar surface with respect to the spin axis of the covere d object 
(for a detailed analysis of the RME see lOhta et aDl2005[) . 

Using a code originally presented in iDreizler et al I (l2Q09l) . 
we have undertaken simulations of the RME for various geomet- 
ric configurations of 2M0535-05 during the primary eclips43 
as it would be seen with the Ultraviolet and Visual Echelle 
Spectrograph (UVES) at the Very Large Telescope (VLT) (see 
Fig-HB- For the data quality we assumed the constraints given 
by the UVES at the VLT exposure time calculatoi0 in version 

3.2.2. The computations show that, using Th-Ar reference spec- 
tra and also the telluric A and B bands as benchmarks, a time 
sampling with one spectrum every 1 245 s and a S/N of > 7 
around 8 600 A are necessary to get 21 measurements during the 
primary eclipse and an accuracy of < 100 m/s. 

In principle, there are four parameters for the background 
object of the transit to be fitted in our simulations of the RME: 
the rotational velocity Vrot, the inclination of the spin axes with 
respect to the line of sight the angle between the projection 
of the spin and the projection of the orbital plane normal onto 
the celestial plane A, and the orbital inclination with respect to 
the line of sight /. From light curve analyses, both rotational ve- 
locities in 2M0535-05 and the orbital inclination / are known. 
Thus, for the simulation of the primary eclipse 7^,2 and X2 are 
the remaining free parameters. 

The obliquities 1 /=i,2, i.e. the real 3 -dimensional angle be- 
tween the orbital normal and the spin axis of the occulted object, 
is related to the other angles as 



cos(i//i) = cos(/^,/) cos(0 -h sin(/^,/) sin(0 cos(Ai). (27) 



^ The 'primary eclipse' refers to the major flux decrease in the sys- 
tem's light curve. Due to the significantly higher effective temperature 
of the secondary mass BD the primary eclipse occurs when the primary 
mass component transits in front of the secondary companion, as seen 
from Earth. 

^ http ://w w w. eso . org/observing/etc 



While the two obliquities are intrinsic angles of the system, 
they cannot be measured directly. They can only be inferred from 
/, i and Ai, which depend on the position of the observer with 
respect to the system. Since we are only interested in the possible 
options for the measurement of th e obliquities in 2M 0535-05, 
we refer the reader to the paper bv lWinnetaD(l2005h for a dis- 
cussion of Eq. (1271) and the geometrical aspects of the RME. 
With / = 88.49° the first term in Eq. (|271) degrades to insignifi- 
cance, which yields cos(</^/) ^ sin(/*,/) cos(Ai). 

At low values for and Ai the fitted solutions to the RME 
are degenerate and there are multiple solutions within a certain 
confidence interval. But our simulations for the transit show that 
the error due to the observational noise is on the same order as 
the error due to degeneracy and thus we find standard deviations 
in Ii, 2 and A2 of cri^ ^ ^ 20° and (T^2 ~ 2^°' respectively. The un- 
certainty in i//2 depends not only on the uncertainties in 1^, 2 and 
A2 but also on the actual values of 7^,2 and A2. But in all cases, 
the standard deviation in the secondary's obliquity (T^2 < 

If present in 2M0535-05, a considerable misalignment of 
the secondary BD of 50° could be detected with a l-cr accuracy 
of 20° or less. Thus, an observed i//2 value of 50° would be a 2.5- 
cr detection of spin-orbit misalignment. Unless RME measure- 
ments suggest i// ^ 90°, RME observations alone are unlikely 
to provide definitive evidence that any of the tidal models we 
consider is responsible for the temperature reversal. 



4.1 .2. Further observations of BD binaries 

Besides the option of RME measurements for testing the geo- 
metric implications, there does exist a possibility to verify our 
estimate of log(2) 3.5 for BDs in general. Comparison of 
observed orbital properties with values constrained by the equa- 
tions that govern the orbital evolution might constrain the free 
parameters, here Q. Using Eq. (O, we find that, assuming only 
a slight initial eccentricity of 0.05, the eccentricity of a BD bi- 
nary system similar to 2M0535-05, in terms of masses, radii, 
rotational frequencies, and semi-major axis would increase to 1 
after ^ 500 Myr if the quality factors of the two BDs are < 10^-^ 
(see left panel in Fig. [12]). A measurement of e in such an evolved 
state could not constrain 2 in a 2M0535-05 analog since either 
the initial eccentricity could have been relatively large while the 
orbit evolved rather slowly due to high Q values or a small ini- 
tial value of e could have developed to a large eccentricity due 
to small values of Q. 

We also simulate the evolution of a 2M0535-05 analog but 
with a diff'erent rotational frequency of the primary constituent in 
order to let the eccentricity decrease with time. We neglected the 
evolution of all the other physical and orbital parameters since 
we are merely interested in a tentative estimate. For a given can- 
didate system the analysis would require a self-consistent cou- 
pled evolution of all the diff'erential equations. For the arbitrary 
case of Pi = P2 = 14.05 d we find that, even for the most ex- 
treme but unrealistic case of an initial eccentricity equal to 1 , this 
fictitious binary would be circularized on a timescale of 100 Myr 
for log(2) < 5 (see right panel in Fig. [12]). Findings of old, ec- 
centric BD binaries with rotational and orbital frequencies that 
yield circularization in the respective system would set lower 
limits to Q. 



4.1.3. Rotational periods in 2l\/l0535-05 

Another, and in fact a crucial, constraint on Q for BDs comes 
from the synchronization time scale 4ynch of the two BDs in 
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2M0535-05. Following the equation given in Lev07 and taking 
the initial orbital mean motion and semi-major axis of the system 
as calculated with an uncoupled system of differential equations 
from model #1 , we derive 4ynch,i =0.07 Myr for the primary and 
4ynch,2 = 0.04 Myr for the secondary with log(2) = 3.5. Since 
the rotation in both BDs is not yet synchronized with the or- 
bit and the age of the system is about 1 Myr, log(2) = 3.5 is not 
consistent with the age of 2M0535-05. Both components should 
have synchronous rotation rates already. We find the consistent 
value for 2 to be > 10^-^, yielding synchronization time scales 
4ynch,i > 0.69 Myr and %nch,2 > 0.37 Myr. 

To make this estimate for Q more robust, we present the evo- 
lution of the BDs' rotational periods in Fig.[T3]and compare it to 
the critical period for a structural breakup Pcrit- The evolution- 
ary tracks are calculated with model #1 and Eq. (30) in FM08. 
As a rough approach we do not couple this equation with those 
for the other orbital parameters. The left panel of Fig. [T3l shows 
that for log(2i) = 3.5 and ij/i = 0° the primary's initial rota- 
tion period 1 Myr ago is ^ 0.3 d. The initial rotation period for 
the secondary, for log(22) = 3.5 and i//2 = 0°, is about -0.2 d, 
where the algebraic sign contributes for a retrograde revolution 
(right panel in Fig. [131). For most of its lifetime, the secondary 
would have had a retrograde rotation and just switched the ro- 
tation direction within the last few 10,000 yr, which is very un- 
likely in statistical terms. Since the orbital momentum is on the 
order of 10^^ kgm^/s and the individual angular momenta are 
about 10"^^ kgm^/s, the shrinking process might not have had a 
serious impact on the rotational evolution. Tides have dominated 
the spin evolu tions. 

Following'S cholz & Eisloffell (l2005l) , the critical breakup pe- 
riod Pent depends only on the body's radius and its mass. The 
radius evolution for BDs is very uncertain for the first Myr after 
formation but we estimate their initial radii to be as large as the 
solar radius. This yields Pcrit,i ~ 0.5 d for both the primary and 
the secondary BD. As stated above, the moduli of the initial ro- 
tation periods of both BDs would have been smaller than 0.5 d 
for Q values of < 10^-^. This inconsistency gives a lower limit 
to Qi and Q2 since values of < 10^-^ would need an initial ro- 
tation periods of both BDs which are smaller than their critical 
breakup periods. Obliquities larger than 0° would accelerate the 
(backwards) evolution and yield even larger lower limits for Qi 
and Q2. Thus, our simulations of the rotational period evolution 
of both BDs require log(2BD) ^ 3.5, whereas the tidal synchro- 
nization timescale even claims log(2BD) ^4.5. 

4.2. Evolutionary embedment of tidal heating 

Tidal heating must be seen in the evolutionary context of 
the system. On the one hand, the tidal energy rates gener- 
ate a temperature increase on the Kelvin-Helmholtz time-scale, 
which is ^ 2 Myr for the BDs in 2M0535-05 - and thus on 
the order of the system's age, as per Eq. ([26b . On the other 
hand, tidal heating will affect the shrinking and cooling pro- 
cess of young BDs in terms of an evo lutionary retardation. As 
models show (ID' Antona & Mazzitellil[i9 97; Bar affe et al. II 19981: 
IChabrier et aril2000l:IChabrier & Baraffdl2000 ). single BDs cool 
and shrink significantly during their first Myrs after formation. 
Adding an energy source comparable to the luminosity of the 
object will slow down the aging processes such that the ob- 
served temperature and luminosity overshoot at some later point 
is not only due to the immediate tidal heating but also due to 
its past evolution. Consequently, the luminosity and temper- 
ature overshoot in the secondary might not (only) be due to 
present-day tidal heating, but it could be a result of an evolu- 



tionary retardation process triggered by the presence of the pri- 
mary as a perturber. Coupled radius-orbit evolutionary models 
have already given plausibl e explanations for the inflated radii 
of some extrasolar planets (iGu et aP l20Q3l: iMiller et al.l l2009l: 
llbgui & Burrowsll2009l:llbgui et al.ll2009allbl) . 

For a consistent description of the orbital and physical his- 
tory of 2M0535-05, one would have to include the evolution of 
obliquities BD radii Ri, eccentricity e, semi-major axis a, and 
rotational frequencies D/. Note that there is a positive feedback 
between radial inflation and tidal heating: as tidal heating in- 
flates the radius, the tidal heating rate can increase and - in turn 
- may cause the radius to inflate even more. In a self-consistent 
orbital and structural simulation of 2M0535-05, tidal inflation, 
neglected in our computations of the Tq^ increase in Eq. ([26b , 
will result naturally from the additional heating term introduced 
by tides. 

In conjunction with 2M0535-05 that means the actual heat- 
ing rates necessary to explain the Tq^ and luminosity excess 
in the secondary are lower than they would have to be if there 
would be no historical context. Relating to Figs.[4l[6l[8l and[TOl 
the implied obliquity and Q factor for the secondary are - again - 
shifted towards lower and higher values, respectively. Embedded 
in the historical context of tidal interaction in 2M0535-05, 
i//2 < 50° and log(22) > 3.5 may also explain the temperature 
reversal and the luminosity excess of the secondary. 

These trends, however, are contrary to that induced by tidal 
inflation. If tidal heating is responsible for a radial expansion 
of 10 and 20% in the primary and secondary, the values of the 
dissipation factor necessary to explain the Teff reversal would be 
^ 0.8 smaller in log(e2) (see Sect.[23b. 

5. Conclusions 

We surveyed four different published tidal models, but neglect 
any evolutionary background of the system's orbits and the com- 
ponents' radii to calculate the tidal heating in 2M0535-05. Our 
calculations based on models #2 and #4, which are most compat- 
ible with the observed properties of the system, require obliqui- 
ties 0, ~ 50° and a quality factor log(2) ^ 3.5 in order 
to explain the luminosity excess of the secondary. Additionally, 
the observed temperature reversal follows naturally since we 
may reproduce a reversal in temperature increase due to tides: 
dT2 > dTi . In model #2, synchronous rotation of the perturbed 
body is assumed. Since this is not given in 2M0535-05, the ac- 
tual heating rates will be even higher than those computed here. 
Our results for the heating rates as per model #2 are thus lower 
limits, which shifts the implied obliquity of the secondary and 
its Q factor to lower and higher values, respectively. 

Considerations of the synchronization time scale for the BD 
duet and the individual rotational breakup periods yield con- 
straints on 2bd for BDs. We derive a lower limit of log(2BD) > 
4.5. This is consistent with estimates of 2- values for M dwarfs, 
log(2dM) ~ 5, and the quality factors of Jupiter, 2 • 10^ < Qq^ < 
2 . 10^ and Neptune, 10^ < log(2g) < lO^-^ (see Sect. [2XTb . 
With log(2BD) > 4.5 tidal heating alone can neither explain the 
temperature reversal in the system nor the luminosity excess of 
the secondary. 

An obliquity of 50°, however, would be reasonable in view 
of recent results from measurements of the RME in several tran- 
siting exoplanet system^ Currently, out of 18 planets there are 
7 with significant spin-orbit misalignments > 30° and some 



^ See http://www.hs.uni-hamburg.de/EN/Ins/Per/Hellerl for an 
overview. 
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of them are even in retrograde orbits around their host stars. 
A substantial obliquity might cause an enhanced heating in 
the 2M0535-05 secondary, while the primary's spin could be 
aligned with the orbital spin, leading to negligible heating in the 
primary. 

Despite the advantages of distance-independent radius and 
luminosity measurements of close, low-mass binaries, the com- 
parison of fundamental properties of the constituents with the- 
oretical models of isolated BDs must be taken with care. This 
applies also to the direct translation from the discrepancies be- 
tween observed and modeled radii for a fixed metallicity into 
an apparent age diff ^erence as a calibration of LMS models 
(IStassun et al.l l2009b . Tidal heating might be a crucial contri- 
bution to discrepancies between predicted and observed radi i 
in other eclipsing low -mass bi nary systems ("Ri bas et aDl2008l) . 
As recently shown by llbgui & Burrows ( 2009), tidal heating in 
extra-solar giant planets in close orbits at a < 0.2 AU with mod- 
est to high eccentricities of ^ > 0.2 can explain the increased 
radii of some planets, when embedded in the orbital history with 
its host star. 

Improvement of tidal theories is necessary to estimate the 
relation between tides and the observed radii of LMS being usu- 
ally too large as compared to models. A tidal model is needed 
for higher orders of arbitrary obliquities and eccent ricities that 
also ac counts for arbitrary rotation rates. As stated by lGreenbergl 
(l2009l) . a formal extension of the simple 'lag-and-add' procedure 
of tidal frequencies the theory of constant phase lag is question- 
able. Besides the extension, conciliation among the various mod- 
els is needed. The results from the models applied here should 
be considered preliminary but are suggestive and indicate the 
possible importance of tides in binary BD systems. 

Several issues remain to be addressed for a more detailed 
assessment of tidal heating in 2M0535-05: /. reconciliation and 
improvement of tidal theories; //. self-consistent simulations of 
the orbital and physical evolution of the system and the BDs; 

measurements of the system's geometric configuration; /v. 
constraints on the tidal quality factors of BDs. 
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Fig. 1. Degeneracy parameter ^ = h^T Kk^T^) (solid line) with 
model parameters similar to those of the 2M0535-05 primary 
and radius-integrated luminosity L (dashed line) as a function of 
radius. To fit into the plot, L is normalized to 10. 
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Fig. 2. Orbital evolution of 2M0535-05 after model #1 going back in time for 1.5 Myr. Left: Eccentricity evolution. Depending 
on Q and on the age of the system, its initial eccentricity has not been smaller than ^ 0.3133, which is ^ 97.4% of its current 
value. Right: Obliquity evolution of the secondary BD for three different values of Q2. Simulations started at 'time = 0' for i//2 ^ 
{0°, 20°, 40°, 60°, 80°} and were evolved backwards in time. For log(22) > 4 there is no significant change in i//2. For all the treated 
values of Q2, the obliquity of the 2M0535-05 secondary is still close to its natal state. 
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Fig. 3. Tidal heating after model #1. Left: (Primary) Projection of E^^^ ^ onto the log(Qi)-i//i plane. The stepsize between contour 
lines is chosen to be A = 0.5 in logC^"*^ ^/W). Right: (Secondary) Projection of E^^^ ^ ^^^^ the log(Q2)-i//2 plane. Although there is a 
dependence on i//2, the energy rates at a fixed value for the quality factor are smaller than those for the primary. 
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Fig. 4. Temperature increase after model #1. Left: (Primary) Projection of dTi onto the log(2i)-j/^i plane. For a significant temper- 
ature increase, Qi would have to be much smaller than 10^-^, but such a temperature increase is not observed in the primary. Right: 
(Secondary) Projection of JT2 onto the log(Q2)-i//2 plane. Even for very low values of Q2 and high obliquities i//2 the observed 
temperature increase cannot be reconstructed. For any given point in the i//-log(Q) plane, dT2 < dTi, which does not support the 
observed temperature reversal. 
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Fig. 5. Tidal heating after model #2. Left: (Primary) Projection of E^^ ^ onto the log(Qi)-i//i plane. The stepsize between contour 
lines is chosen to be A = 0.5 in log(£^*^ ^/W). The tidal energy rates strongly depend on a putative obliquity, different from model 
#1. Right: (Secondary) Projection of E^^ ^ onto the log(22)-<A2 plane. For the three models (#1, #2, and #4) invoking Q and i//, these 
rates are the highest of all for any given point in the i/^-\og(Q) plane - for the primary as well as for the secondary. 
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Fig. 6. Temperature increase after model #2. Left: (Primary) Projection of dTi onto the log(2i)-0^i plane. Right: (Secondary) 
Projection of JT2 onto the log(22)-<A2 plane. For any given location in the \og(Q)-i// plane, model #2 yields the strongest tem- 
perature increase compared to the other models - both for the primary and the secondary, respectively. For a given spot in Q-i/r 
space there is an inversion in temperature increase: dr2 > dTi, i.e. the less massive BD is heated more. 
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Fig. 7. Tidal heating within the primary (soHd Hne) and secondary (dashed Hne) after model #3. While the tidal heating rate of the 
secondary becomes comparable to its observed luminosity overshoot for r > 100 s, if the same r is applied to the primary, heating 
within the primary would lead to a larger luminosity than is observed. 
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Fig. 8. Temperature increase of the primary (solid line) and secondary (dashed line) after model #3. Contrary to what is observed, 
the primary would be hotter than the secondary. 
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Fig. 9. Tidal heating after model #4. Left: (Primary) Projection of E^^ ^ onto the log(Qi)-i//i plane. The stepsize between contour 
lines is chosen to be A = 0.5 in log(E^^^/W). Right: (Secondary) Projection of E^^^ onto the log(Q2)-i//2 plane. 
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Fig. 10. Temperature increase after model #4. Left: (Primary) Projection of dTi onto the \og(Qi)-i//i plane. Right: (Secondary) 
Projection of JT2 onto the log(22)-<A2 plane. For log(22) ~ 3.5 and an obliquity of i//2 ~ 70° the temperature increase of the 
secondary becomes similar to the observed one. For the whole range of Q and i// there is an inversion in temperature increase, 
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Primary eclipse duration ~ 7h:31min (±50min) Primary eclipse duration ~ 7h:31min (± 50min) 



Fig. 11. Simulations for the Rossiter-McLaughlin effect as it would be seen with UVES during the primary eclipse of 2M0535-05, 
which occurs when the secondary mass BD is occulted by the primary. The S/N is 7. Left: The orbital inclination / is fixed at 88.49° 
(see Table [B and /I = 0, which means the transiting primary BD follows a path parallel to the secondary's equator. The alignment 
of the secondary's spin axis 7^,2 varies between 90° (perpendicular to the line of sight) and 50°. Right: With / fixed at 88.49° and 
h,2 = 90°, A2 varies between 0° (primary path parallel to the secondary's equator) and 60° (primary path strongly misaligned with 
the secondary's equator). 
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Fig. 12. Orbital evolution of a 2M0535-05 analog after model #1. Left: Eccentricity evolution for different values of Q for the next 
500 Myr. The initial eccentricity was arbitrarily chosen: eini = 0.05. For log(2) ^3.5 this binary will be disrupted within 500 Myr. 
Right: Eccentricity evolution of a 2M0535-05 analog but with Pi = P2 = 14.05 d for different values of Q. Contrary to the scenario 
in the left figure, the changed rotational period of the primary BD now leads to circularization of the system. Measurements of e in 
LMS binaries with known ages can give lower limits to Q. 
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Fig. 13. Rotational evolution of the two BDs in 2M0535-05 after model #1 for different values of Qi and Q2. Left: (Primary) Going 
backwards in time, the rotation period decreases. For log(2i) = 3.5, Pi drops below the critical period for structural breakup of 
^ 0.5 d already before the date of birth around 1 Myr ago. Right: (Secondary) For log(22) = 5.5 we show the tracks for i//2 = 0° 
and 80° for comparison. For log(22) = 4.5 the rotation direction switches at about -0.1 8 Myr and for log(22) = 3.5 at roughly 
-10, 000 yr. 



